Multi-dimensional pattern recognition by use of digital holography

ABSTRACT

A method and system for performing three-dimensional pattern recognition by use of digital holography is disclosed. The complex amplitude distribution generated by an object is recorded by phase-shifting interferometry. The digital hologram contains information about the objects shape, location, and orientation. This information allows one to perform multi-dimensional pattern-recognition methods with a high degree of discrimination and to measure three-dimensional orientation changes.

TECHNICAL FIELD

[0001] This invention relates multi-dimensional object, pattern or data recognition by use of digital holography.

BACKGROUND OF THE INVENTION

[0002] Optical correlation methods have proved to be very useful in the design of two-dimensional (2D) pattern-recognition applications (A. VanderLugt, IEEE Trans. Inf. Theory IT-10, 139 (1964); J. L. Homer and P. D. Gianino, Appl. Opt. 23, 812 (1984); D. Psaltis, E. G. Paek, and S. S. Venkatesh, Opt. Eng. 23, 698 (1984); Ph. Refrefgier, Opt. Lett. 15, 854 (1990) and J. W. Goodman, Introduction to Fourier Optics (McGraw-Hill, New York, 1996) which are incorporated herein by reference in their entirety).

[0003] Recently, there has been increasing interest in three-dimensional (3D) optical information processing because of its vast potential applications. Several methods have been proposed to extend optical correlation methods to three-dimensional object recognition. One approach consists of processing different 2D projections of both a three-dimensional scene and a three-dimensional reference object by use of conventional 2D Fourier methods (A. Pu, R. Denkewalter, and D. Psaltis, Opt. Eng. 36, 2737 (1997) which is incorporated herein by reference in its entirety).

[0004] Other methods also involve the acquisition of different 2D perspectives, but the recognition is performed by three-dimensional Fourier-transform methods (J. Rosen, Opt. Lett. 22, 964 (1997) which is incorporated herein by reference in its entirety). This approach has analogies to incoherent-light three-dimensional imaging based on spatial coherent functions (H. Arimoto, K. Yoshimori, and K. Itoh, Opt. Commun. 170, 319 (1999) which is incorporated herein by reference in its entirety). Fourier-transform profilometry has also been applied to three-dimensional shape recognition (J. J. Esteve-Taboada, D. Mas, and J. Garcia, Appl. Opt. 22, 4760 (1999) which is incorporated herein by reference in its entirety).

[0005] However, holography, seems to be a more attractive method of performing three-dimensional image recognition, since a single hologram is able to record three-dimensional information on the object, avoiding sequential recording of several 2D perspectives and maintaining the phase (H. J. Caulfield, ed., Handbook of Optical Holography (Academic, London, 1979) which is incorporated herein by reference in its entirety).

SUMMARY OF THE INVENTION

[0006] A method and system for recognizing multi-dimensional objects, patterns or sets of data utilizing digital holography is disclosed. The method comprises generating a hologram of a reference object; generating a hologram of an input object; and correlating the hologram of the reference object with the hologram of the input object generating thereby a set of correlation peaks. The presence of a correlation peak indicates similarity between the reference object and the input object and the lack of the presence of a correlation peak indicates dissimilarity between the reference object and the input object. The method may be used to determine a change in a kinematic property of an object.

[0007] A storage medium is disclosed that is encoded with a set of data created by generating a hologram of a reference object, generating a hologram of an input object, and correlating the hologram of the reference object with the hologram of the input object generating thereby a set of correlation peaks.

[0008] The system for performing pattern recognition comprises a light source for generating an optical beam; means for dividing the optical beam into a reference beam and an object beam; means for introducing a phase shift between the reference beam and the object beam; an object positioned within the object beam; a beam combiner for combining the reference beam and the object beam; and a detector for detecting the combination of the reference beam and the object beam.

EXPLANATION OF THE DRAWINGS

[0009]FIG. 1 is a schematic representation of a phase-shifting interferometer for recording an inline digital hologram of a three dimensional object.

[0010]FIG. 2 is a schematic representation of the relationship between the displacement of a hologram window and the angle of view of the hologram.

[0011]FIG. 3A is a first reconstruction of a three dimensional object from a digital hologram in a plane orthogonal to a detector.

[0012]FIG. 3B is a second reconstruction of a three dimensional object from a digital hologram in a plane orthogonal to a detector.

[0013]FIG. 4A is the autocorrelation by digital holography of the reconstructed three dimensional object represented in FIG. 3A.

[0014]FIG. 4B is the cross correlation by digital holography of the three-dimensional objects represented in FIGS. 3A and 3B.

[0015]FIG. 5A is a conventional two-dimensional autocorrelation of the image of the two dimensional object in FIG. 3A.

[0016]FIG. 5B is a conventional two-dimensional cross correlation of the images of the two-dimensional objects in FIGS. 3A and 3B.

[0017]FIG. 6A is the correlation peak versus the angle of view used to measure three dimensional object orientation for the correlation of the three-dimensional object in FIG. 3A with a rotated version of the same three-dimensional object obtained by digital holography.

[0018]FIG. 6B is the correlation peak versus the angle of view for the correlation of the three-dimensional object of FIG. 3A with a rotated version of the same three-dimensional object giving the maximum correlation peak.

[0019]FIG. 7 is a schematic representation of a computer network connected to the phase-shifting holographic system of FIG. 1.

[0020]FIG. 8 is a representation of the four step method of inducing a phase shift in the reference beam of the interferometric system of FIG. 1.

[0021]FIG. 9A is a first schematic representation of an optical correlation of a hologram of a reference object and a hologram of an input object.

[0022]FIG. 9B is a second schematic representation of an optical correlation of a hologram of a reference object and a hologram of an input object.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0023] An optical system for the digital holographic recording of a three-dimensional object, a pattern or a set of data is based on an interferometer, such as a Mach-Zehnder interferometer 100, and is depicted in FIG. 1. Polarized light 104 emanating from a laser 102, such as an Argon laser, a helium-neon laser an infrared laser or other such laser, is divided into an object beam 104 b and a reference beam 104 a. The object beam 104 b illuminates a still or moving three-dimensional reference object 130, while the reference beam 104 a forms an in-line interference pattern together with the light 104 d diffracted by the reference object 130. The reference object 130 may comprise, for example, an optical image, a digitized image, a computer generated image, a one-dimensional or multi-dimensional object, pattern or set of data, an electrical signal or an optical signal. The object 130 may also include two or three-dimensional phase objects or color objects and is shown as a three-dimensional object by way of exemplification.

[0024] The interference pattern, or hologram, is formed or recorded on a detector or storage medium 128, such as a CCD camera. In FIG. 1, the reference numeral 108 denotes a spatial filter, 110 denotes a lens, 112 is a diaphragm, 114 and 126 are beam splitters, 106, 120, 122, 124, are plane mirrors, 116 and 118 are phase retarders with retardation of λ/2 and λ/4 respectively, where λ is the wavelength of the light emanating from the laser 102. The reference object 130 is located at a distance d from a recording plane 210 in the CCD camera 128. The CCD camera 128 is in turn in signal communication with a monitor 136 and a computer 134.

[0025] Consider the opaque three-dimensional reference object 130, as a complex amplitude distribution U_(O)(x, y, z), where x, y are the transverse coordinates and z denotes the paraxial distance from the reference object 130 to the output plane 210 of the CCD camera 128. In the Fresnel approximation, neglecting secondary diffraction, the complex amplitude distribution, H_(O)(x, y,), due to the object beam 104 d and recorded in the output plane 210 can be evaluated from the following superposition integral: $\begin{matrix} \begin{matrix} {{H_{O}\left( {x,y} \right)} = \quad {{A_{H}\left( {x,y} \right)}{\exp \left\lbrack {{\varphi}_{H}\left( {x,y} \right)} \right\rbrack}}} \\ {= \quad {\frac{1}{i\quad \lambda}{\int{\int{\int_{\infty}^{\infty}{{U_{O}\left( {x^{\prime},{y^{\prime};z}} \right)}\frac{1}{z}{\exp \left( {2\frac{\pi}{\lambda}\quad z} \right)} \times}}}}}} \\ {\quad {{\exp \left\lbrack {{\frac{\pi}{\lambda \quad z}\left\lbrack {\left( {x - x^{\prime}} \right)^{2} + \left( {y - y^{\prime}} \right)^{2}} \right\rbrack}} \right\rbrack}{x^{\prime}}{y^{\prime}}{z}}} \end{matrix} & (1) \end{matrix}$

[0026] The complex field generated by the reference beam 104 a at the output plane 210 can be written as R(x, y; α_(n))=A_(R)exp[i(φ+α_(n))], where φ is the constant phase when both fast axes of the phase retarders 116, 118 are aligned with the direction of polarization and α_(n) is the phase shift in other configurations of the phase retarders 116, 118. As seen in FIG. 8, by adjusting the relative position of the slow and fast axes of the phase retarders 116, 118, phase shifts of α_(n)=0, π/2, π, and 3π/2 are produced without altering the constant amplitude, A_(R), of the reference beam 104 a. The phase retarders 116, 118 allow for the modulation or shifting of the phase, α_(n), of the reference beam 104 a with respect to the object beam 104 b. In this manner, phase shifts, α_(n), are introduced into the reference beam 104 a with respect to the object beam 104 b by the relative positioning of the fast and slow axes of the retarder plates 116, 118 with the polarization of the incident light 104 a. The aforesaid phase shifts, α_(n), may also be introduced into the reference beam 104 a with respect to the object beam 104 b by other methods or apparatus, such as the use of a liquid crystal chip in place of the phase retarders 116, 118 or by the movement of mirror 124.

[0027] The complex amplitude given by Eq. (1) is measured by recording four interferometric patterns, or holograms, I_(n)(x, y; α_(n))=|H_(O)(x, y)+R(x, y; α_(n))|².

[0028] The holograms, I_(n)(x, y; α_(n)), formed by the combination of H_(O)(x, y) and R(x, y; α_(n)) are recorded digitally, or on film, as a hologram at the output plane 210. Such recording may be in the form, for example, of tangible storage media, such as optical storage materials (e.g., photorefractive materials), or digitally on floppy diskettes, CD-ROMs, hard drives, optically or electrically addressable spatial light modulators, charge coupled devices or any other computer-readable storage medium or processing unit addressable across a distributed environment such as a computer network system. It will be recognized that the holograms formed at the output plane 210 are not limited to the in-line holograms of FIG. 1, but may also be formed by off-axis holograms as is well known in the art.

[0029] The phase difference between the object beam 104 b and the reference beam 104 a is given by $\begin{matrix} {{{\varphi_{H}\left( {x,y} \right)} - \phi} = {\arctan \left\lbrack \frac{{I\left( {x,{y;{{- 3}{\pi/2}}}} \right)} - {I\left( {x,{y;{{- \pi}/2}}} \right)}}{{I\left( {x,{y;0}} \right)} - {I\left( {x,{y;{- \pi}}} \right)}} \right\rbrack}} & (2) \end{matrix}$

[0030] The amplitude of the interference pattern of the object beam 104 d and the reference beam 104 a can be calculated from: $\begin{matrix} {{{A_{H}\left( {x,y} \right)}A_{R}} = {{\frac{1}{4}{I\left( {x,{y;0}} \right)}} - \frac{I\left( {x,{y;{- \pi}}} \right)}{\left. \left. {\cos\left\lbrack {{\varphi_{H}\left( {x,y} \right)} - \phi} \right.} \right) \right\rbrack}}} & (3) \end{matrix}$

[0031] where the cosine function is obtained from Eq. (2). The parameters φ and A_(R) are constant factors which can be replaced with 0 and 1, respectively.

[0032] Different regions or windows within the digital hologram of the reference object 130, record different perspectives, views or segments of the reference object 130. These different perspectives, views or segments of the reference object can be reconstructed by the numerical propagation of H_(O)(x, y). The discrete complex amplitude distribution of the reconstructed reference object, U_(O)′(m, n; a_(x), a_(y)), at any plane in the object beam 104 d, within the paraxial approximation and within the angle limitations imposed by the size of the CCD 128, can be numerically computed, aside from constant factors, by the discrete Fresnel transformation: $\begin{matrix} {{U_{O}^{\prime}\left( {m,{n;a_{x}},a_{y}} \right)} = {{\exp \left\lbrack {\frac{\pi}{\lambda \quad d^{\prime}}\left( {{\Delta \quad x^{2}m^{2}} + {\Delta \quad y^{2}n^{2}}} \right)} \right\rbrack} \times {\sum\limits_{m^{\prime} = 0}^{N - 1}\quad {\sum\limits_{n^{\prime} = 0}^{N - 1}\quad {{H_{O}^{\prime}\left( {m^{\prime},{n^{\prime};a_{x}},a_{y}} \right)}{\exp \left\lbrack {\frac{\pi}{\lambda \quad d^{\prime}}\left( {{\Delta \quad x^{\prime 2}m^{\prime 2}} + {\Delta \quad y^{\prime 2}n^{\prime 2}}} \right)} \right\rbrack}{\exp \left\lbrack {- {{2\pi}\left( {\frac{m\quad {m!}}{N_{x}} + \frac{{nn}!}{N_{y}}} \right)}} \right\rbrack}}}}}} & (4) \end{matrix}$

[0033] where d′=−d; (m′, n′) are discrete spatial coordinates in the plane of the hologram 210; (m, n) correspond to discrete spatial coordinates in the plane of the reconstructed object and N_(x) and N_(y) are the number of samples in the x and y directions in the plane of the reconstructed object, respectively. The spatial resolution in the plane of the hologram 210 is (Δx′, Δy′), and the spatial resolution in the plane of the reconstructed object is (Δx, Δy). H_(O)′(m, n, a_(x), a_(y)) denotes the amplitude distribution of the object over a window 212 defined within the hologram 210 and used for the reconstruction of the reference object 130 (FIG. 2). H_(O)′(m, n; a_(x), a_(y)) is modified by a linear phase factor, exp└i2π(a_(x)m+a_(y)n)┘to take into account the direction of observation. Thus, $\begin{matrix} {{H_{O}^{\prime}\left( {m,{n;a_{x}},a_{y}} \right)} = {{H_{O}\left( {m,n} \right)}{{rect}\left( {\frac{m - a_{x}}{b_{x}},\frac{n - a_{y}}{b_{y}}} \right)} \times {\exp \left\lbrack {{2\pi}\left( {{a_{x}m} + {a_{y}n}} \right)} \right\rbrack}}} & (5) \end{matrix}$

[0034] wherein rect(g, h) is the rectangle function. In Eqs. (4) and (5), (a_(x), a_(y)) are the pixel coordinates of the window 212 and (b_(x),b_(y)) denote the transverse size of the window 212. FIG. 2 shows the relationship between the displacement of the hologram window 212 and the angle of view along the direction of observation 204. Angles of γ and β are equal to a_(x)Δx′/d and a_(y)Δy′/d, respectively. The minimum angle, γ_(min), between different perspectives of the reference object 130 depends on the CCD pixel size, Δx′, and the distance d, such that γ_(min)=Δx′/d. A similar relationship applies for the y direction for the angle β, i.e., β_(min)=Δx′/d. The maximum angle, γ_(max), is a function of the CCD pixel size and the window size, such that γ_(max,)=(N_(x)−b_(x))Δ_(x)′/d; and similarly for the angle β in the y direction, i.e., β_(max,)=(N_(y)−b_(y))Δy′/d. It will be appreciated that the resolution of the reconstructed images, Δx and Δy are also affected by the window size, i.e., Δx=λd/(b_(x)Δx′) and Δy=λd/(b_(y)Δy′). Equation (4) can be numerically evaluated with a fast Fourier-transform method. An alternative method of evaluating Eq. (4) is to use the propagation transfer function method: $\begin{matrix} {{U_{O}^{\prime}\left( {m,{n;a_{x}},a_{y}} \right)} = {F^{- 1}\left( {{F\left\lbrack {H_{O}\left( {m,{n;a_{x}},a_{y}} \right)} \right\rbrack} \times {\exp \left\lbrack {{\pi\lambda}\quad {d^{\prime}\left\lbrack {\frac{u^{2}}{\left( {\Delta \quad x\quad N_{x}} \right)^{2}} + \frac{v^{2}}{\left( {\Delta \quad y\quad N_{y}} \right)^{2}}} \right\rbrack}} \right\rbrack}} \right)}} & (6) \end{matrix}$

[0035] where F denotes the fast Fourier transformation and u and v are discrete spatial frequencies.

[0036] To perform recognition of, or discrimination between, three-dimensional objects, patterns or sets of data, a Fourier-matched filter approach is applied to the information obtained by digital holography. Consider a three-dimensional input object 132 whose complex amplitude distribution is given by U_(P)(x, y, z). The input object 132 may comprise an optical image, a digitized image, a one dimensional set of data, a two dimensional set of data, a multi-dimensional set of data, an electrical signal or an optical signal and is shown as a three-dimensional object by way of exemplification. The input object 132 is also located at a distance d from the CCD 128, and has a real term in a corresponding Fresnel hologram 210, of H_(P)(x, y). Correlation between different views of the input object 132 and a given view of the reference object 130, i.e., object recognition, can be performed by defining a window, H_(P)′(x, y; a_(x)′, a_(y)′), within the Fresnel hologram, H_(P)(x, y), of the input object 132 by the use of Eq. (5), followed by reconstructing the input object 132 with Eq. (4) (yielding the reconstructed input object U_(P)′(x, y, z)) and computing the correlation between U_(P)′(x, y, z) and U_(O)′(x, y, z). Alternatively, by use of Eq. (6), we can write the correlation intensity of the reference amplitude, U_(O)′(x, y; a_(x), a_(y)), with that of the input, U_(P)′(x, y; a_(x)′, a_(y)′), generated from the holograms: $\begin{matrix} \begin{matrix} {{C_{OP}\left( {x,{y;a_{x}},{a_{y};a_{x}^{\prime}},a_{y}^{\prime}} \right)} = \quad {{F^{- 1}\left\{ {{F\left\lbrack {U_{O}^{\prime}\left( {x,{y;a_{x};a_{y}}} \right)} \right\rbrack} \times} \right.}}} \\ \left( {{\quad \left. {F*\left\lbrack {U_{P}^{\prime}\left( {x,{y;a_{x}^{\prime}},a_{y}^{\prime}} \right)} \right\rbrack} \right\}}} \right)^{2} \\ {= \quad {F^{- 1}\left\{ {{F\left\lbrack {H_{O}^{\prime}\left( {x,{y;a_{x}},a_{y}} \right)} \right\rbrack} \times} \right.}} \\ \left( {{\quad \left. {F*\left\lbrack {H_{P}^{\prime}\left( {x,{y;a_{x}^{\prime}},a_{y}^{\prime}} \right)} \right\rbrack} \right\}}} \right)^{2} \end{matrix} & (7) \end{matrix}$

[0037] Thus, by performing the correlation between different regions or windows defined within the hologram of the reference object 130 and the hologram of the input object 132, properly modified by a linear phase factor, an evaluation of the correlation between different perspectives, views or segments of the reference object 130 and input object 132 is performed. It will be appreciated that rough objects involve fast fluctuations of the reconstructed phase during translations, thus reducing the shift invariance. Nevertheless, high sensitivity is obtained for small object motion. This is illustrated by measuring a small rotation of a rough three-dimensional object.

[0038] A three-dimensional object recognition experiment was performed with two reproductions of cars with an approximate size of 25 mm×45 mm. The objects were located at an approximate distance, d, of about 865 mm from the output plane 210. The pictures in FIGS. 3A and 3B were obtained by recording a digital hologram of each car, denoted H_(O)(x, y) and H_(P)(x, y). Reconstruction was performed by computing Eq. (4) with a fast Fourier-transform method applied to both of the holograms.

[0039]FIG. 4A shows a correlation peak 230 in a plot 232 a of the autocorrelation of the reconstructed object shown in FIG. 3A, which was performed by use of a digital hologram, H_(O)′(x, y; 0, 0), obtained by placement of a window with size b_(x)=b_(y)=256 pixels centered at the origin of H_(O)(x, y) (see Eq. (5)). FIG. 4B shows a plot 232b of the cross correlation of the three-dimensional object in FIG. 3A with that in FIG. 3B, which was also obtained by digital holography. No correlation peak is seen in FIG. 4B. The reference hologram is the previous digital hologram H_(O)′(x, y; 0, 0), and the input, H_(P)′(x, y; 0, 0), is obtained by placement of a window with the same size centered at the origin of H_(P)(x, y). Both plots are normalized to the same value. By way of comparison, FIG. 5A shows the conventional two-dimensional autocorrelation 232 c and FIG. 5B the conventional two-dimensional cross correlation 232 d, using the objects in FIGS. 3A and 3B registered as two-dimensional intensity images.

[0040] In FIG. 9A the correlation of the hologram 222 of the reference object 130 with the hologram 224 of the input object 132 may be performed optically rather than digitally. An electrically addressable spatial light modulator 220 a is interfaced with the detector 128. The digital holograms 222, 224 of the reference and input objects are displayed on the electrically addressable spatial light modulator 220 a and the correlation of the holograms is constructed in a correlation plane by a reference wave 104 g. In FIG. 9B, the holograms 222, 224 are formed by an optically addressable spatial light modulator 220 b such as a liquid crystal light valve or a ferro-electric spatial light modulator. The hologram 224 of the input object 132 and the hologram 222 of the reference object 130 are read out from the output of the optically addressable spatial light modulator 220 b by the reference wave 104 g and the correlation of the holograms 222, 224 is formed at the correlation plane 232. The following are incorporated herein by reference: B. Javidi and J. L. Horner, “Real-time Optical Information Processing,” Academic Press, 1994; U.S. Pat. No. 4,832,447, entitled Joint Transform Image Correlation Using a Nonlinear Spatial Light Modulator at the Fourier Plane; U.S. Pat. No. 4. 5,841,907, entitled Spatial Integrating Optical Correlator for Verifying the Authenticity of a Person, Product or Thing; U.S. Pat. No. 5,699,449, entitled Method and Apparatus for Implementation of Neural Networks for Face Recognition; U.S. Pat. No. 5,367,579, entitled Method of Removing Spurious Responses from Optical Joint Transform Correlators; U.S. Pat. No. 5,119,443, entitled Nonlinear Joint Transform Optical Correlator Having Varying Degrees of Nonlinearity; U.S. Pat. No. 5,485,312, entitled Optical Pattern Recognition System and Method for Verifying the Authenticity of a Person, Product or Thing; B. Javidi, “Optical Information Processing,” in Encyclopedia of Electrical and Electronic Engineering, Volume on Optics, John Wiley, 1999 and B. Javidi, “Optical Neural Networks,” in Encyclopedia of Electrical and Electronic Engineering, Volume on Neural Networks, John Wiley, 1998 and B. Javidi and E. Tajahuerce, “Three Dimensional Image Processing And Reconstruction” International Symposium On Photonics For Aerospace Application Of Optics, SPIE Proceedings Vol. 4043, Orlando Fla., Apr. 24 - 28, 2000

[0041] It will be appreciated that the methods disclosed herein can be used to measure changes in the kinematic properties (translation or rotation) of the three-dimensional input object 132 with respect to the reference object 130 by adjusting the parameters in the phase factor applied to the window, H_(P)′(x, y; a_(x)′,a_(y)′), defined within the Fresnel hologram, H_(P)(x, y), of the input object 132. In particular, FIG. 6 shows the correlation obtained by digital holography between the three-dimensional object in FIG. 3A, which was used as a reference, and the same slightly rotated three-dimensional object used as input. A window with a size, for example of, b_(x)=b_(y)=256 pixels and centered at the origin of the hologram of the reference object was utilized. The same window is shifted to different positions in the hologram of the input object. FIG. 6A shows the maximum value of the correlation 232 e as a function of the displacement (a_(x),a_(y)) of the window in the hologram of the input. A peak 232 e is obtained for a_(x)=−12, a_(y)=−2, corresponding to an angle of view (γ,β)=(0.007°, 0.001°). FIG. 6B shows the correlation 232 f for this perspective of the three-dimensional object. The sharp peak 232 f indicates that the object is only rotated with respect to the axial view of the reference.

[0042] Thus a method of determining a change in a kinematic property of an object, is realized by generating a first hologram of the object, generating a second hologram of the object, correlating the first hologram with the second hologram, generating thereby a set of correlation peaks and analyzing the set correlation peaks. The aforesaid method may be accomplished by defining a window within the first hologram, defining a window within the second hologram and correlating the first and second windows, generating thereby a set of correlation peaks and analyzing the set correlation peaks. The aforesaid method may further comprise autocorrelating the first hologram and comparing the autocorrelation of the first hologram with the correlation of the first hologram with the second hologram.

[0043] Referring to FIG. 7, the detector 128 is shown connected to a computer or communications network 400, such as a local area network (LAN) or a wide area network (WAN) via a computer 134. The detector 128 may also be connected directly to a liquid crystal display (LCD), liquid crystal television (LCTV) or an electrically or optically addressable spatial light modulator or ferroelectric spatial light modulator 412 a. The computer network 400 includes a plurality of similar client personal computers 404 connected to a server 410 from remote geographical locations by wired or wireless connections, by radio based communications, by telephony based communications, or by other network-based communications. The computer 134 may also be connected directly to another like computer 134 or to a display device 412, such as a liquid crystal display (LCD), liquid crystal television (LCTV) or an electrically or optically addressable spatial light modulator (SLM) for 3D optical reconstruction of the objects 130, 132 or patterns or sets of data. The computer network 400 is in turn similarly connected to other computers 502 or networks 504 through the Internet 500. The computers 134, 404, 502, display devices 136, 412, 412 a and other devices of the networks 400, 500, 504 are configured to execute program software, that allows them to send, receive, store, record and process the holograms, H_(O)(x, y) and H_(P)(x, y), either separately or in combination, between and amongst themselves via the networks 400, 504 or the Internet 500. Such processing includes, for example, image compression and decompression, filtering, contrast enhancement, image sharpening, noise removal and correlation for image classification. Decompressed images may be displayed on display devices such as liquid crystal displays, liquid crystal TVs or electrically or optically addressable spatial light modulators.

[0044] Based upon the foregoing description, an optoelectronic holographic method and system for performing pattern recognition of three-dimensional objects, or sets of data, by use of digital phase-shifting holography has been presented. This method is based on the ability of a digital plane hologram to reconstruct different perspectives of a three-dimensional object. The method allows one to obtain three-dimensional information on reference and the input objects. Three-dimensional pattern recognition is carried out by use of a digital matched filter technique applied to the holographic information. The holographic data can be transmitted through conventional digital communications channels to remote locations and reconstructed therefrom digitally or optically.

[0045] Utilizing phase-shifting interferometry, (J. H. Bruning, D. R. Herriott, J. E. Gallagher, D. P. Rosenfeld, A. D. White, and D. J. Brangaccio, Appl. Opt. 13, 2693 (1974) and I. Yamaguchi and T. Zhang, Opt. Lett. 22, 1268 (1997) which are incorporated herein by reference in their entirety) the methods herein electronically record the complex amplitude distribution generated by a three-dimensional object at a single plane located in the Fresnel diffraction region. Similarly, a digital hologram of a three-dimensional reference pattern is recorded to be used as a correlation filter. Thus, pattern recognition by use of three-dimensional information can be performed by application of correlation methods to digital holograms.

[0046] While preferred embodiments have been shown and described, various modifications and substitutions may be made thereto without departing from the spirit and scope of the invention. Accordingly, it is to be understood that the present invention has been described by way of illustration only, and such illustrations and embodiments as have been disclosed herein by reference are not to be construed as limiting the claims. 

What is claimed is:
 1. A method of pattern recognition comprising: generating a hologram of a reference object; generating a hologram of an input object; and correlating the hologram of the reference object with the hologram of the input object to generate a set of correlation peaks.
 2. The method as set forth in claim 1 further comprising analyzing the set of correlation peaks generated by the correlation of the hologram of the reference object with the hologram of the input object wherein the presence of a correlation peak indicates similarity between the reference object and the input object and wherein the lack of the presence of a correlation peak indicates dissimilarity between the reference object and the input object.
 3. The method as set forth in claim 1 further comprising recording the hologram of the reference object.
 4. The method as set forth in claim 1 further comprising recording the hologram of the input object.
 5. The method as set forth in claim 1 wherein the correlating of the hologram of the reference object with the hologram of the input object comprises: transforming the hologram of the reference object; transforming the hologram of the input object; and multiplying the transformation of the hologram of the reference object by the transformation of the hologram of the input object.
 6. The method as set forth in claim 5 wherein the transforming of the hologram of the reference object comprises Fourier transforming the hologram of the reference object and wherein transforming the hologram of the input object comprises Fourier transforming the hologram of the input object.
 7. The method as set forth in claim 6 further comprising transforming the multiplication of the transformation of the hologram of the reference object and the transformation of the hologram of the input object.
 8. A method of pattern recognition comprising: generating a hologram of a reference object; generating a hologram of an input object; defining a window within the hologram of the input object; correlating the hologram of the reference object with the window defined within the hologram of the input object to generate a set of correlation peaks.
 9. The method as set forth in claim 8 further comprising analyzing the set of correlation peaks generated by the correlation of the hologram of the reference object with the window defined within the hologram of the input object wherein the presence of a correlation peak indicates similarity between the reference object and the input object and wherein the lack of the presence of a correlation peak indicates dissimilarity between the reference object and the input object.
 10. The method as set forth in claim 8 further comprising recording the hologram of the reference object.
 11. The method as set forth in claim 8 further comprising recording the hologram of the input object.
 12. The method as set forth in claim 8 wherein the correlating of the hologram of the reference object with the window within the hologram of the input object comprises: transforming the hologram of the reference object; transforming the window defined within the hologram of the input object; and multiplying the transformation of the hologram of the reference object by the transformation of the window defined within the hologram of the input object.
 13. The method as set forth in claim 12 wherein the transforming of the hologram of the reference object comprises Fourier transforming the hologram of the reference object and wherein transforming the window within the hologram of the input object comprises Fourier transforming the window within the hologram of the input object.
 14. The method as set forth in claim 13 further comprising transforming the multiplication of the transformation of the hologram of the reference object and the transformation of the window defined within the hologram of the input object.
 15. The method as set forth in claim 8 wherein the defining of a window within the hologram of the input object comprises defining a window according to the equation ${{H_{p}^{\prime}\left( {m,{n;a_{x}},a_{y}} \right)} = {{H_{p}\left( {m,n} \right)}{{rect}\left( {\frac{m - a_{x}}{b_{x}},\frac{n - a_{y}}{b_{y}}} \right)}}},$

wherein H_(P)′(m, n; a_(x),a_(y)) is the defined window, ${rect}\left( {\frac{m - a_{x}}{b_{x}},\frac{n - a_{y}}{b_{y}}} \right)$

is the rectangle function, H_(P)(m, n) is the hologram of the reference object, a_(x) is the location of the center of the window in the x direction, a_(y) is the location of the center of the window in they direction, b_(x) is width of the window in the x direction, and b_(y) is width of the window in the y direction.
 16. A storage medium encoded with a set of data created by: generating a hologram of a reference object; generating a hologram of an input object; and correlating the hologram of the reference object with the hologram of the input object to generate a set of correlation peaks.
 17. The storage medium as set forth in claim 16 further comprising analyzing the set of correlation peaks generated by the correlation of the hologram of the reference object with the hologram of the input object wherein the presence of a correlation peak indicates similarity between the reference object and the input object and wherein the lack of the presence of a correlation peak indicates dissimilarity between the reference object and the input object.
 18. The storage medium as set forth in claim 16 wherein the generating of the hologram of the reference object includes recording the hologram of the reference object.
 19. The storage medium as set forth in claim 16 wherein the generating of the hologram of the input object includes recording the hologram of the input object.
 20. The storage medium as set forth in claim 16 wherein the correlating of the hologram of the reference object with the hologram of the input object comprises: transforming the hologram of the reference object; transforming the hologram of the input object; and multiplying the transformation of the hologram of the reference object by the transformation of the hologram of the input object.
 21. The storage medium as set forth in claim 20 wherein the transforming of the hologram of the reference object comprises Fourier transforming the hologram of the reference object.
 22. The storage medium as set forth in claim 20 wherein the transforming of the hologram of the input object comprises Fourier transforming the hologram of the input object.
 23. A storage medium encoded with a set of data created by: generating a hologram of a reference object; generating a hologram of an input object; defining a window within the hologram of the input object; correlating the hologram of the reference object with the window within the hologram of the input object to generate a set of correlation peaks.
 24. The storage medium as set forth in claim 23 further comprising analyzing the set of correlation peaks generated by the correlation of the hologram of the reference object with the hologram of the input object wherein the presence of a correlation peak indicates similarity between the reference object and the input object and wherein the lack of the presence of a correlation peak indicates dissimilarity between the reference object and the input object.
 25. The storage medium as set forth in claim 23 wherein the generating of the hologram of the reference object includes recording the hologram of the reference object.
 26. The storage medium as set forth in claim 23 wherein the generating of the hologram of the input object includes recording the hologram of the input object.
 27. The storage medium as set forth in claim 23 wherein the correlating of the hologram of the reference object with the window within the hologram of the input object comprises: transforming the hologram of the reference object; transforming the window within the hologram of the input object; and multiplying the transformation of the hologram of the reference object by the transformation of the window within the hologram of the input object.
 28. The storage medium as set forth in claim 27 wherein the transforming of the hologram of the reference object comprises Fourier transforming the hologram of the reference object.
 29. The storage medium as set forth in claim 27 wherein the transforming of the window within the hologram of the input object comprises Fourier transforming the window within the hologram of the input object.
 30. The storage medium as set forth in claim 23 wherein the defining of a window within the hologram of the input object comprises defining a window according to the equation ${{H_{p}^{\prime}\left( {m,{n;a_{x}},a_{y}} \right)} = {{H_{p}\left( {m,n} \right)}{{rect}\left( {\frac{m - a_{x}}{b_{x}},\frac{n - a_{y}}{b_{y}}} \right)}}},$

wherein H_(P)′(m, n; a_(x),a_(y)) is the defined window, ${rect}\left( {\frac{m - a_{x}}{b_{x}},\frac{n - a_{y}}{b_{y}}} \right)$

is the rectangle function, H_(P)(m, n) is the hologram of the reference object, a_(x) is the location of the center of the window in the x direction, a_(y) is the location of the center of the window in the y direction, b_(k) is width of the window in the x direction, b_(y) is width of the window in the y direction.
 31. The method as set forth in claim 7 further comprising taking the absolute value of the transformation of the multiplication of the transformation of the hologram of the reference object by the transformation of the hologram of the input object.
 32. The method as set forth in claim 31 further comprising squaring the absolute value of the transformation of the multiplication of the transformation of the hologram of the reference object by the transformation of the hologram of the input object according to the equation C _(OP)(x, y; a _(x) , a _(y) ; a _(x) ′, a _(y)′)=|F ⁻¹ {F[H _(O)′(x, y; a _(x) , a _(y))]×F*[H _(P)′(x, y; a _(x) ′, a _(y)′)]}|² wherein C_(OP)(x, y; a_(x), a_(y); a_(x)′, a_(y)′) is the intensity of the correlation of the hologram of the reference object with the hologram of the input object, H_(O)′(x, y; a_(x), a_(y)) is the hologram of the reference object, H_(P)′(x, y; a_(x)′, a_(y)′) is the window defined within the hologram of the input object, F indicates a transformation, F¹ indicates an inverse transformation and F* indicates a complex transformation.
 33. The method as set forth in claim 14 further comprising taking the absolute value of the transformation of the multiplication of the transformation of the hologram of the reference object by the transformation of the hologram of the input object.
 34. The method as set forth in claim 33 further comprising squaring the absolute value of the transformation of the multiplication of the transformation of the hologram of the reference object by the transformation of the hologram of the input object according to the equation C _(OP)(x, y; a _(x) , a _(y) ; a _(x) ′, a _(y)′)=|F ⁻¹ {F[H _(O)′(x, y; a _(x) , a _(y))]×F*[H _(P)′(x, y; a _(x) ′, a _(y)′)]}|² wherein C_(OP)(x, y; a_(x), a_(y); a_(x)′, a_(y)′) is the intensity of the correlation of the hologram of the reference object with the hologram of the input object, H_(O)′(x, y; a_(x), a_(y)) is the hologram of the reference object, H_(P)′(x, y; a_(x)′, a_(y)′) is the window defined within the hologram of the input object F indicates a transformation, F¹ indicates an inverse transformation and F* indicates a complex transformation.
 35. The method as set forth in claim 32 wherein F indicates a Fourier transformation.
 36. The method as set forth in claim 34 wherein F indicates a Fourier transformation.
 37. The method as set forth in claim 3 wherein the recording of the hologram of the reference object includes digitally recording the hologram of the reference object.
 38. The method as set forth in claim 37 wherein the digitally recording of the hologram of the reference object comprises storing the hologram of the reference object in a computer readable storage medium.
 39. The method as set forth in claim 4 wherein the recording of the hologram of the input object includes digitally recording the hologram of the input object.
 40. The method as set forth in claim 39 wherein the digitally recording of the hologram of the input object comprises storing the hologram of the input object in a computer readable storage medium.
 41. The method as set forth in claim 10 wherein the recording of the hologram of the reference object includes digitally recording the hologram of the reference object.
 42. The method as set forth in claim 41 wherein the digitally recording of the hologram of the reference object includes storing the hologram of the reference object in a computer readable storage medium.
 43. The method as set forth in claim 11 wherein the recording of the hologram of the input object includes digitally recording the hologram of the input object.
 44. The method as set forth in claim 43 wherein the digitally recording of the hologram of the input object includes storing the hologram of the input object in a computer readable storage medium.
 45. The method as set forth in claim 15 further comprising applying a phase factor, exp└i2π(a_(x)m+a_(y)n)┐, to the window defined within the hologram of the input object.
 46. A method of determining a change in a kinematic property of an object, the method comprising: generating a first hologram of the object; generating a second hologram of the object; and correlating the first hologram with the second hologram, to generate a set of correlation peaks.
 47. The method as set forth in claim 46 further comprising analyzing the set of correlation peaks wherein the presence of a correlation peak indicates a change in a kinematic property of the object and wherein the lack of the presence of a correlation peak indicates no change in the change in a kinematic property of the object.
 48. A method of determining the change in a kinematic property of an object, the method comprising: generating a first hologram of the object; defining a first window within the first hologram of the object; generating a second hologram of the object; defining a second window within the second hologram of the object; correlating the first window and the second window to generate a set of correlation peaks.
 49. The method as set forth in claim 46 further comprising autocorrelating the first hologram and comparing the autocorrelation of the first hologram with the correlation of the first hologram with the second hologram.
 50. The method as set forth in claim 48 further comprising analyzing the correlation peaks wherein the presence of a correlation peak indicates a change in a kinematic property of the object and wherein the lack of the presence of a correlation peak indicates no change in the kinematic property of the object.
 51. The method as set forth in claim 1 wherein the reference object comprises an optical image, a digitized image, a one dimensional set of data, a two dimensional set of data, a multi-dimensional set of data, an electrical signal, an optical signal, a two-dimensional phase object, a multi-dimensional phase object or a color object.
 52. The method as set forth in claim 1 wherein the input object comprises an optical image, a digitized image, a one dimensional set of data, a two dimensional set of data, a multi-dimensional set of data, an electrical signal, an optical signal, a two-dimensional phase object, a multi-dimensional phase object or a color object.
 53. The method as set forth in claim 1 further comprising conveying the hologram of the reference object to a remote location over a distributed computer network.
 54. The method as set forth in claim 1 further comprising conveying the hologram of the input object to a remote location over a distributed computer network.
 55. The method as set forth in claim 1 further comprising conveying the correlation of the hologram of the reference object with the hologram of the input object to a remote location over a distributed computer network.
 56. The method as set forth in claim 8 further comprising conveying the hologram of the input object to a remote location over a distributed computer network.
 57. The method as set forth in claim 8 further comprising conveying the hologram of the reference object to a remote location over a distributed computer network.
 58. The storage medium as set forth in claim 8 further comprising conveying the correlation of the hologram of the reference object with the window defined within the hologram of the input object to a remote location over a distributed computer network.
 59. The storage medium as set forth in claim 23 further comprising conveying the hologram of the input object to a remote location over a distributed computer network.
 60. The storage medium as set forth in claim 23 further comprising conveying the hologram of the reference object to a remote location over a distributed computer network.
 61. The storage medium as set forth in claim 23 further comprising conveying the correlation of the hologram of the reference object with the window within the hologram of the input object to a remote location over a distributed computer network.
 62. The method as set forth in claim 1 further comprising processing the hologram of the reference object by image compression of the hologram.
 63. The method as set forth in claim 62 further comprising conveying the compressed hologram to remote locations over a distributed computer network.
 64. The method as set forth in claim 1 further comprising processing the hologram of the input object by image compression of the hologram.
 65. The method as set forth in claim 64 further comprising conveying the compressed hologram to remote locations over a distributed computer network.
 66. The method as set forth in claim 62 further comprising processing the compressed hologram of the reference object by image decompression of the hologram.
 67. The method as set forth in claim 64 further comprising processing the compressed hologram of the input object by image decompression of the hologram.
 68. A system for performing pattern recognition, the apparatus comprising: a light source for generating an optical beam; means for dividing the optical beam into a reference beam and an object beam; means for introducing a phase shift between the reference beam and the object beam; an object positioned within the object beam; a beam combiner for combining the reference beam and the object beam; and a detector for detecting the combination of the reference beam and the object beam.
 69. The system as set forth in claim 68 wherein the detector comprises a computer addressable detector.
 70. The system as set forth in claim 69 wherein the detector is connected to a distributed computer network.
 71. The system as set forth in claim 70 wherein the detector is connected to an electrically or optically addressable spatial light modulator.
 72. The system as set forth in claim 68 wherein the object comprises a two-dimensional or three-dimensional phase object, a color object, an original set of data comprising an optical image, a digitized image, a computer generated image, a one dimensional set of data or multi-dimensional set of data, an electrical signal or an optical signal.
 73. The system as set forth in claim 68 wherein the light source is a source of infrared light.
 74. The method as set forth in claim 8 wherein the reference object comprises an optical image, a digitized image, a one dimensional set of data, a two dimensional set of data, a multi-dimensional set of data, an electrical signal, an optical signal, a two-dimensional phase object, a multi-dimensional phase object or a color object.
 75. The method as set forth in claim 8 wherein the input object comprises an optical image, a digitized image, a one dimensional set of data, a two dimensional set of data, a multi-dimensional set of data, an electrical signal, an optical signal, a two-dimensional phase object, a multi-dimensional phase object or a color object.
 76. The method as set forth in claim 46 wherein the reference object comprises an optical image, a digitized image, a one dimensional set of data, a two dimensional set of data, a multi-dimensional set of data, an electrical signal, an optical signal, a two-dimensional phase object, a multi-dimensional phase object or a color object.
 77. The method as set forth in claim 46 wherein the input object comprises an optical image, a digitized image, a one dimensional set of data, a two dimensional set of data, a multi-dimensional set of data, an electrical signal, an optical signal, a two-dimensional phase object, a multi-dimensional phase object or a color object.
 78. The method as set forth in claim 48 wherein the reference object comprises an optical image, a digitized image, a one dimensional set of data, a two dimensional set of data, a multi-dimensional set of data, an electrical signal, an optical signal, a two-dimensional phase object, a multi-dimensional phase object or a color object.
 79. The method as set forth in claim 48 wherein the input object comprises an optical image, a digitized image, a one dimensional set of data, a two dimensional set of data, a multi-dimensional set of data, an electrical signal, an optical signal, a two-dimensional phase object, a multi-dimensional phase object or a color object.
 80. The system as set forth in claim 68 wherein the detector comprises a photorefractive material.
 81. The method as set forth in claim 1 wherein correlating the hologram of the reference object with the hologram of the input object comprises: displaying the holograms of the input object and the reference object on electrically or optically addressable spatial light modulators; addressing the electrically or optically addressable spatial light modulators with a reference beam; and processing the holograms optically to generate a set of correlation peaks.
 82. The method as set forth in claim 8 wherein correlating the hologram of the reference object with the window defined within the hologram of the input object comprises: displaying the window defined within the hologram of the input object and the hologram of the reference object on electrically or optically addressable spatial light modulators; addressing the electrically addressable spatial light modulator with a reference beam; and processing the holograms optically to generate a set of correlation peaks.
 83. The method as set forth in claim 1 wherein correlating the hologram of the reference object with the hologram of the input object comprises: forming the holograms of the input object and the reference object on electrically or optically addressable spatial light modulators; addressing the optically addressable spatial light modulator with a reference beam; and processing the holograms optically to generate a set of correlation peaks.
 84. The method as set forth in claim 8 wherein correlating the hologram of the reference object with the window defined within the hologram of the input object comprises: forming the window defined within the hologram of the input object and the hologram of the reference object on electrically or optically addressable spatial light modulators; addressing the electrically or optically addressable spatial light modulators with a reference beam; and processing the holograms optically to generate a set of correlation peaks. 